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Conventionally, one calculates a zero in a beta function by computing this function to a given loop 
order and solving for the zero. Here we discuss a different method which is applicable in theories 
where one can perform a partial diagrammatic summation to infinite-loop order. We show that this 
method, compared with the conventional method, yields much better agreement with exact results 
in the case of an asymptotically free gauge theory with Af = 1 supersymmetry. Applications to 
other field theories are also discussed. 


In a quantum field theory, the dependence of the in¬ 
teraction coupling on the Euclidean momentum scale, /i, 
where it is measured is of fundamental importance. This 
dependence is described by the beta function, /?, of the 
theory [Ij. Of particular interest is a zero of P, since at 
this zero, the coupling is scale-invariant. We focus on a 
zero of p at nonvanishing coupling. The (physical) zero 
of P closest to the origin in coupling constant space is 
an infrared (IR) zero for an asymptotically free theory 
and an ultraviolet (UV) zero for an IR-free theory. In 
a conventional analysis of the beta function of a given 
theory, one considers the series expansion of /? to a given 
order in powers of the coupling and investigates whether 
or not it has a (physical) zero away from the origin in 
coupling-constant space. 

Here we discuss a different approach to the calculation 
of a zero of the beta function. Our approach is appli¬ 
cable to theories where one can perform a summation, 
to infinite-loop order, of the terms from a (sub)set of 
diagrams contributing to the beta function. Although 
our method can be applied to any theory where one can 
carry out this type of summation, it is conveniently illus¬ 
trated in the context of an asymptotically free vectorial 
gauge theory. Thus, we consider such a gauge theory, 
with gauge group G, running coupling g = g(p), and Nf 
massless Dirac fermions in a representation R of G Q. 
Let a(/i) = g(/i) 2 /(47t). If Nf is sufficiently large (still 
maintaining asymptotic freedom), then the two-loop beta 
function has an IR zero. As g decreases from large val¬ 
ues in the UV, a(fi) increases toward the value of this 
zero, which thus plays an important role in determin¬ 
ing the IR properties of the theory Q. The beta func¬ 
tion is p g = dg/dt, where dt = d\ng , or equivalently, 
Pa = da/dt = [g/(2n)]/3 g . This has the series expansion 


p a = - 2 a ^2 b e a 1 , ( 1 ) 

l=i 

where he, is the Uloop coefficient and a = g 2 /(16n 2 ). The 
coefficients b\ and 62 are independent of the regulariza¬ 
tion scheme, while the be for £ > 3 are scheme-dependent 
f4|. The property of asymptotic freedom means bi > 0. 
The n-loop ( n£ ) beta function, denoted P a ,ni, is given by 
Eq. (JT]) with i = n as the upper limit on the sum. We 


define the reduced (r) function 


P 


a,r — 


Pa 

[a 2 b 1 /(2n)\ 



( 2 ) 


and Pa, r ,nii given by Eq. fl2J) with the upper limit t = n 
on the sum. The conventional approach to the study of 
a possible IR zero in P a , n e consists of an analysis of the 
n — 1 zeros of the n-loop reduced beta function p a ,r,nt- 
In particular, at the two-loop level, the IR zero is given 
by aiR t 2 i = —47 t&i/& 2, which is physical if 62 < 0. Since 
we are interested in an IR zero here, we thus assume that 
the matter content of the theory is such that 62 < 0. 

However, if one can carry out a summation of an infi¬ 
nite (sub)set of diagrams contributing to the beta func¬ 
tion, this yields a result that, although formally equiva¬ 
lent to m when expanded in a series, involves a closed- 
form (c/) functional factor, / c /. Normalizing this factor 
so that f c f = 1 at a = 0, we can reexpress p a in this case 
as 


Pa = 



( 3 ) 


i.e., Pa,r = fcf fs, where f s is given by a series (s) ex- 
pansion 


fs = l + '£fs, j ai . (4) 

3 =1 


In a theory where this summation can be carried out 
and p a can thus be expressed in the form ([3]), our pro¬ 
posed method is to analyze the zeros in f c f and f s rather 
than the zeros of P a ,r,ni as a means of gaining informa¬ 
tion about the physical zero in the beta function. This 
method takes advantage of the information from the all¬ 
orders summation, while the conventional method of an¬ 
alyzing zeros in P a ,r,ni to a given finite loop order n does 
not. In particular, if f c f has no physical zero, then the 
proposed method is to examine the zeros in f s . 

Let us consider an illustrative function for f c f, namely 


* cf 1 + dia ’ 


( 5 ) 


(where d\ is a constant), so that P a , r = fs/{ 1 + did). 
Then 


fs — 1 + -g- ^2 fin + ^1 bn- 1 ) a 


n —1 


n=2 


( 6 ) 
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Solving the equation f s = 0 to the lowest nontrivial order, 
n = 2, we obtain for the IR zero 

47T 

= -(bM + d, ’ (7) 

which is physical if ( 62 /& 1 ) + di < 0 (the subscript cfs 
connotes the direct dependence on f s and the implicit 
dependence on f c f-) Since we assume b 2 < 0 (and b\ > 
0 ), a-iR'Cfs is physical if 


appeal that a number of exact results are known for it 0 ]- 
11 ], enabling one to compare approximate calculations 
with these exact results 0 ] - 0 - We will show that the 
method discussed here, compared with the conventional 
method, yields much better agreement with exact results 
and solves a puzzle found in [131 ] . 

In this SGT the constraint of asymptotic freedom re¬ 
quires that Nr < Nf^iz, where Nf^iz = 3Ca/( 2T/) [Il- 
it is known 0-0 that if Nf is in the interval Inacp 
defined by Nf tCr < Nf < Nf^u, where 


di 



( 8 ) 


Nf,ar 


3 C A 
AT f 


Nf,blz 
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( 12 ) 


Because the existence of the IR zero in f) a is scheme- 
independent, one may require that the scheme used 
for the summation to infinite-loop order that yields ([3]) 
should not remove this IR zero, and hence that d\ satis¬ 
fies the inequality (0). As compared with aiR >2 e, 


_ OlIR^t 

&IR,cfs _ ^ aiR,2i " 

1 4-7T 


(9) 


Hence, if d\ is negative (positive), then ajR iC f s is smaller 
(larger) than oiir^i- As is evident from the dependence 
on d\ , aiR^cfs incorporates information from the sum¬ 
mation to infinite-loop order, in contrast to ajR tU e for 
any finite n. A necessary condition for the result ([9]) to 
be physically meaningful is that if f c f has any divergent 
physical singularities, they must occur farther from the 
origin than In our example, f c f has a simple pole 

at ctcf^poie = — 4n/di, which occurs at physical coupling 
if d\ < 0. The required inequality is a c / iPO z e > aiR tC f s , 
i.e., 

11 

d\ (b 2 /bi) + di 


Since &2 < 0 for the existence of olir^h and we only need 
to address the case where d\ is negative, (IH is equivalent 
to the inequality 


1 1 
W\ > T\wM+W\' 


( 11 ) 


which is clearly true. This proves that (for the case we 
consider here, where the theory has an IR zero in fd a , 2 e) 
even if f c f contains a physical pole, one has a c / jPO z e > 
aiR^cfs, so this pole is irrelevant for the evolution from 
weak coupling in the UV to the IR zero at aj_R )C /s- 
To demonstrate the usefulness of the proposed method, 
we apply it to an asymptotically free vectorial Af = 1 su¬ 
persymmetric gauge theory (denoted SGT) with gauge 
group SU(iV c ) and Nf pairs of massless chiral super- 
fields Qj, * 1 >j, 1 < j < Nf transforming according to 
the representations R and R of SU(IV C ), respectively. In 
this theory, the scheme-independent coefficients in f3 a are 
h = 3 C A - 2 T f N f and b 2 = 6 C\ - 4 (C A + 2 C f )T f N f 
s i- Calculations to three-loop Q order have been 
done, e.g., in the DR scheme 0. This theory has the 


then this theory flows from weak coupling in the UV 
to a conformal non-Abelian Couloumb phase (NACP) 
in the IR without any spontantaneous chiral symmetry 
breaking (S\SB). A quantity of interest is the anomalous 
dimension of the fermion bilinear, jl 6 j . This has the 
series expansion 


7 m = ^2 Ct - al > 
e =1 


(13) 


where c\ = ACf is scheme-independent, and the q with 
t > 2 are scheme-dependent. Using instanton methods, 
Novikov, Shifman, Vainshtein, and Zakharov (NSVZ) ob¬ 
tained an exact closed-form solution for f} a 0,0 (in a 
particular scheme 0 ): 


/3a,NSVZ 


oF_ 61 - 2T f N nm 
2tt [ 1 - 2C'a a 


(14) 


As Nf decreases from Nf biz & n d approaches Nf cr from 
above, increases from 0 to the maximum value allowed 
by unitarity in this conformal phase, namely 7 m = 1 
jla . 0 . Applying the condition that should saturate 
this unitarity upper bound at an IR zero of /? as Nf 
approaches Nf^ cr from above is one way to derive the 
value of Nf tCr . 

Using known results for the be and eg for 1 < £ < 3 
1-0, Ref. 0] calculated aiR :U e at the n = 2 and 
n = 3 loop level and evaluated the n-loop anomalous 
dimension of the fermion bilinear at a = aiR <n g, denoted 
7 iR,nt- Although higher-order calculations of aiR i7l g are 
scheme-dependent, they are quite valuable in obtaining 
more accurate information about 7 m at an exact IR fixed 
point, which is a universal quantity. Ref. 0] found that 
at the two-loop level, 7 /^, 2 ^ increases monotonically as 
Nf decreases from Nf^iz, but exceeds unity at a value 
of Nf e Inacp larger than Nf tCr . As was noted in 0, 
this is understandable in view of the fact that the two- 
loop perturbative calculations of (3 and 7 m are not ex¬ 
pected to agree precisely with exact results; the differ¬ 
ences provide a quantitative measure of the accuracy of 
these two-loop calculations. Normally, one would expect 
that calculating a quantity to higher-loop order should 
give a more accurate result if perturbative computations 
are reliable. However, the study of the IR behavior at 
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three-loop level in El yielded a different and very puz¬ 
zling result: 'fiR^t does not approach more closely to 
1 as Nf decreases through Inacp toward Nf tCr \ on the 
contrary, 7 /^ 3 ^ reaches a maximum at a small positive 
value and then decreases, passing through zero to nega¬ 
tive values. 

We next show how our proposed method removes this 
puzzling behavior and produces excellent agreement with 
exact results for this SGT. First, observe that p a ,NSVZ 
has the form of Eq. d3j) with 


and 


fcf,NSVZ 


1 

1 - 2 C A a ’ 


fs,NSVZ = 1 — 


2TfNfJm 

bi 


(15) 


(16) 


We focus on the range of Nf where P a ,NSVZ has an IR 
zero. Applying our method, we observe that f c f,NSVZ 
has no zero, so we solve the equation f s ,NSVZ = 0. Re¬ 
taining the maximal scheme-independent information in 
7 m , namely the i = 1 term in m, this equation reads 
b\ — 8NfTfCfa = 0, whence 


_ tt(3C a - 2 T f Nf) 
aiR, c f s 2 T f CfNf 


Cf 


'2 Nf >cr 


- 1 


(17) 


As Nf decreases from Nf^iz to Nf tCr , aiR, c f s increases 
monotonically from 0 to a c f s = n/Cf. Since d\ < 0 
(and with b- 2 < 0 so that otmpi is physical), it follows, 
as a special case of ((9]), that, am tC fs < &ir,u in this 
SGT. Our general discussion above shows that the pole 
in f c f,NSVZ at a r f r,ni* nsvz = 27 t/Ca lies farther from 
the origin than aiR iC f s [20]. As Nf approaches Nf b i~ 
from below, both aiR tC f s and chir^i approach zero, and 
the ratio aiR, c f s /oLiR, 2 t /* 1 . 

Thus, aiR'Cfs retains perturbative reliability much 
better than the conventional two-loop result 


am,2i 


2tt(3CU - 2TfNf) 

2(C A + 2C f )T f N f -3C* A ’ 


(18) 


which diverges as Nf \ Nf^ b2z = 3C\/[2(C A + 2C/)T/], 
where b 2 —> 0. For R equal to the fundamental represen¬ 
tation (with Young tableau □), Nf tb2z G Inacp- 

Next, we study y m , calculated to its maximal scheme- 
independent order, = 4(7/a, and evaluated at 
a iR,cfs- We obtain 


to be a monotonically increasing function with decreasing 

Nf € Inacp El- 

Second, 7 /fl iC / s < 1 for all Nf £ Inacp , in agreement 
with the constraint from unitarity. In contrast, 'fiRfii = 
was found to violate this unitarity constraint for a range 
of Nf toward the lower end of the interval Inacp [131 ]: 
indeed, it was found that jiR^e diverges, just as uiRfu 
does, as Nf \ Nf :b2z [l3[. As was noted in 13], this 
behavior indicated that the perturbative calculations of 
o-ir. 21 . and the resultant evaluation of 7 m at this value of 
c-ir .21 are not reliable in this region. With our different 
procedure, we avoid this pathological behavior. 

Third, we find that 


7 iR,cfs /* 1 as Nf \ Nf iCr . (20) 


Thus, this 7 /_r, c /s has precisely the three properties of the 
exact all-orders calculation of 7 ir. Therefore, we have 
shown that our method of calculating the IR zero of the 
beta function as a zero of f s rather than as the zero of 
Pa,r,ne, together with the evaluation of 7 m calculated to 
its maximal scheme-independent order to obtain 7 iR, c fs, 
yields excellent agreement with exact results and avoids 
the problems with 7 /^, 2 ^ and 7 /^, 3 ^ found in fl3| . 

For the case R = O Eqs. (ED and ED are 


_ 2ttN c (3N c - N f ) 
a iR,cfsD - (7V c 2 _ l)N f 

and 

_ 31V C 

'Tm.cfsD ~ N f 1 ' 


( 21 ) 


( 22 ) 


In this R = □ case, it is also of interest to consider the 
’t Hooft-Veneziano limit N c —> 00 and Nf -A 00 with 
r = Nf/N c fixed and finite, and a(/j,)N c a finite function 
of /i. In this limit, Inacp is3/2<r<3. In terms of 
£ = aN c in this limit, we have 

^ 2tt(3 — r) 3 -r 

^IR,cfs ~ - , 7 IR,cfs ~ - ■ (23j 


In contrast, 


I,IR,2l 


27t( 3 — r) 
2r — 3 


(24) 


which diverges as r \ 3/2 at the lower boundary of 
Inacp , and 


JIR,cfs 


3CU _ IV/,& iz _ . _ 2Nf tCr 

2 TfNf N f N f 


1 . (19) 


")IR,2t 


r(i — 1)(3 — r) 
2(3 - 2r ) 2 ’ 


(25) 


From ED. it follows that 7 m \ 0 as Nf Nf :b i z - Fur¬ 
thermore, from (11911 we obtain three important results 
for 7 m , which exhibit complete agreement with exact 
results and avoid the problems with 'yiR^e and 
found in [r|. First, Eq. (fTiil) shows that as Nf decreases 
from Nf b i z to Nf cr in the interval Inacpi hR,cfs in¬ 
creases monotonically. In contrast, 7 /_r, 3 ^ was found not 


which violates its unitarity upper bound as r decreases 
through r = 2 in the interior of Inacp [HD E3 and di¬ 
verges as r \ 3/2. 

Since the bf, with l > 3 are scheme-dependent, so is 
the closed-form function f c f resulting from the infinite- 
loop summation. One might think that the analysis of an 
IR zero in (3 a could be simplified by applying a scheme 
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transformation (ST) to eliminate the be with £ > 3. An 
explicit ST that reduces j3 a to the two-loop expression 
in the local neighborhood of the origin, a = 0 , was con¬ 
structed and studied in [2lj]-f23j (see also IH). However, 
it was shown that STs that are acceptable for small a of¬ 
ten lead to unphysical results when applied at a generic 
zero of a beta function away from the origin 21|-[23l|. 
We find this to be true for STs that attempt to remove 
be with l > 3 in the present SGT. Hence, one must deal 
with the full series CD. where the method discussed here 
greatly improves the analysis. 

Another demonstration of the usefulness of the method 
discussed here is obtained from analysis of a test function 
for fd a . A basic question for a study of the zero(s) of a 
beta function and, more generally, in mathematics is: 
how well do the zero(s) of a truncation of a power series 
for a function f(z) reproduce the zero(s) of the exact 
function. In Ref. [221] this question was investigated for 
a general asymptotically free gauge theory, using for j3 a>r 
the test function 


Pa,r = h{a) = 


sin(7rvd ) 


(7rvd ) 


(26) 


where a = a/am. We analyzed the zeros of truncations 
of the Taylor series expansion of h(a ) to high order and 
determined how rapidly these approach am, i.e., how 
rapidly the zeros in the ratios a approach 1. Here, to 
apply our method, we start from the above series for 
h(a) and perform an infinite-order summation of part of 
it, using the identity (sin a:)/# = [( 2 /x) sin(x/ 2 )] cos(x/ 2 ) 
with x = 7 ryfdt, taking ( 2 /x) sin(a:/ 2 ) = f c f to be the 
closed-form function and taking cos(a:/ 2 ) = /„, expanded 
as a series. That is, we write 


sin a; 
x 


y, (—l ) 71 X 2n 

h ( 2 -+ 1 )! 


sin(a;/ 2 ) 

x/2 


yv (-1)" {X/2) 2U 
^ ( 2 nj! 


(27) 

This f c f has no zero in the interval 0 < x < n, i.e., 
0 < a < am, so the IR zero in /5 a ^ r comes from f s . 
We calculate this from the truncation of the series for f s 
and denote it as a c / s>n . In Table Q] we list the values of 
S n = am^ n — 1 obtained by (i) the conventional method 
of analyzing I3 a ^ r ,n as a series, as in [ 22 ], denoted {5 n )h, 
and by (ii) the method discussed here, with ( 5 n ) c fs = 
0-cf s,n — 1- As is evident, our present method yields faster 
convergence toward the exact result. 

For a non-supersymmetric gauge theory (in d = 4 di¬ 
mensions) with Nf fermions in a general representation R 
of the gauge group G, a rigorous closed-form expression 
for /3 a has not, so far, been obtained, although a conjec¬ 
ture for an all-orders beta function has been presented in 
[25l . Higher-loop analyses of an IR zero in /3 a have been 
performed using the conventional method of studying the 
zeros of the polynomial /3 a ,r,ni up to the four-loop level 
[26l l27i . These obtained the expected result that calcu¬ 
lating am tn £ and 'yiR^i to higher-loop order led to more 
stable results. For example, for G = SU(3) and Nf = 12 


TABLE I: Values of ( 8 n )h and ( 8n) c f s as a function of respective 
series truncation order n. 


n 

(8n)h 

( Sn)cfs 

2 

-0.392 

-0.189 

3 

complex 

2.78 x 10“ 2 

4 

-3.97 x 10“ 2 

-1.13 x 10 -3 

5 

4.52 x 10 -3 

3.15 x 10 -5 

6 

-2.83 x 10" 4 

-0.592 x 10 -6 

7 

1.35 x 10 -5 

0.805 x 10” 8 

8 

-0.493 x 10~ 6 

-0.829 x 10 _1 ° 


fermions in R = Q = 0.77, ^ir,m = 0.31 and 

1ir,U = 0.25 ( 27 )]. The four=loop result agrees well with 
the lattice measurements - 7 = 0.27(3) and 7 m ~ 0.25 
[13 and 7 = 0.235(46) @0- 

There are also IR-free theories where one can obtain 
an exact closed-form result for the beta function with 
a UV zero at nonzero coupling. A first example is the 
N — > 00 limit of the nonlinear O (N) a model in d = 2 + e 
dimensions 31]. Here, for small e, one finds the closed- 
form result [311] 

*-£-'*-;)■ (28 > 

where A is the effective coupling, A c = 2?re/N; x = 
lirn/v-^oo AIV, and x c = 2i re is the UV zero of f) x . In 
the notation of Eq. ©, here f c f = 1 — ( x/x c ) and 
f s = 1. A second example is supersymmetric quantum 
electrodynamics (SQED) with Nf pairs of chiral super¬ 
fields with charges that can be taken to be 1 and — 1 . 
As is well known, (3 a for this theory can be obtained 
from (fl4l) by setting Ca = 0, Tf = 1, and Cf = 1, 
and is ]3j| Pc,sqed = [a 2 Nf/(2n)\ [1 + (a/br)]. Up to 
the reversal of sign, this has the form of Eq. m with 
fcf = 1 + (a/n) and and f s = 1. It has no UV zero. 
Third, one may consider the SGT with Nf > Nf t bi z . 
Using (ED again for this theory, one observes that the 
numerator of ED is — 1 61 1 — 2NfTf"/ m . Evaluating 7 m 
at the maximal scheme-independent level, 7 m = 4C/a, 
it follows that this numerator has no UV zero. Some re¬ 
lated work for a nonsupersymmetric U(l) gauge theory 
(in d = 4) is in [33j, [3J| ; in this theory, both a conven¬ 
tional analysis of up to the five-loop level and use 
of some partial closed-form results also led to evidence 
against a UV zero [Hj . 

In summary, we have discussed a method for calculat¬ 
ing a zero in the beta function of a quantum field the¬ 
ory that is applicable in cases where one can perform an 
infinite-loop summation of a subset of diagrams. We have 
shown that this method, compared with the conventional 
fixed-loop-order method, yields much better agreement 
with exact results in the case of an asymptotically free 
gauge theory with Af = 1 supersymmetry. We have also 
commented on the application of this method to some 
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IR-free theories. 
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